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1. Introduction 

To understand the physics of how magnetic flux enters a type-II superconductor under 
the influence of apphed currents and flelds is important not only for fundamental reasons 
but also for potential applications. It is well known that magnetic flux enters in the form 
of quantized vortices [1], each carrying magnetic flux 0o = h/2e concentrated within 
an area of the order of A^, where A is the London penetration depth. In the case of 
greatest practical interest, when the intervortex spacing is much smaller than A, the 
macroscopic magnetic flux density B (averaged over a few intervortex spacings) has 
magnitude B = n(/)o, where n is the local areal density of vortices. 

The orientation of B serves as a local reference direction. The Lorentz force per 
unit volume on the vortex array is F = J x B. When F exceeds the local pinning force 
density F^, the vortices generally move with a velocity v in the direction of F^ thereby 
producing an electric fleld [2] E — B x v. In most cases of practical interest, such as 
in superconducting magnets, the magnetic induction in the superconductor produced 
by the self-flelds is perpendicular to the current direction. To calculate quantities of 
practical interest, such as ac losses, one may then use the well-known critical-state 
theory [3, 4, 5, 6, 7], which assumes that B and J are perpendicular to each other 
and E is perpendicular to B. This critical-state theory characterizes the threshold for 
depinning in terms of a critical current density Jc± = Fp/B^ where here we use the 
subscript _L as a reminder that this refers to the case when J is perpendicular to B. 

It has been known for many years, however, that when a type-II superconducting 
wire in a parallel magnetic fleld carries a current above the wire's critical current, a 
flnite electric fleld appears along the length of the wire [8]. Inside the wire, there must 
be a place where J, E^ and B are parallel, not perpendicular. This parallel component 
of the electric fleld cannot be understood in terms of the collective motion of an array 
of unbroken vortices but can be understood as a consequence of flux-line cutting, the 
intersection and cross-joining of locally nonpar allel vortices [5, 8]. This is the analogue 
of the process of vortex-line reconnection in turbulent superfluid ^He, recently fllmed 
by Paoletti et al [9]. Flux-line cutting can be triggered by the helical vortex expansion 
instability [10], in which initially straight vortex lines subjected to a sufflciently large 
current density are susceptible to the growth of helical perturbations. In this paper we 
use Jc\\ to denote the critical current at the threshold for flux-line cutting, where we use 
the subscript || as a reminder that this refers to the case when J is parallel to B. 

Although the behaviour of type-II superconductors is well understood for the case 
that J is perpendicular to the theory for the case that J has components both 
perpendicular and parallel to B is still in a state of development. Such a theory is 
needed, for example, to calculate the ac losses in power transmission cables consisting 
of hehcally wound layers of 2G coated-conductor composite tapes [11]. If the tapes 
are tightly wrapped with no gaps between them, the self-fleld generated inside the 
cable has azimuthal and longitudinal components but no radial components. In such 
a geometry, the vortices in the tapes are subjected to a local current density with not 
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only perpendicular components, which can drive the vortices either inward or outward 
in the radial direction, but also parallel components, which can lead to local helical 
instabilities and subsequent flux-line cutting [10]. It is the purpose of this paper to 
contribute towards the development of a critical-state theory accounting for both flux 
transport and flux cutting. 

In Sees. 2, 3, 4, 5, and 6 we discuss the different predictions of five theoretical 
critical-state models for how the critical-current density Jc and the electric field E (for 
J just above Jc) should vary with the angle between J and B. While there has been 
previous experimental work on fiux cutting in thin films with applied magnetic fields 
tilted out of plane but without in-plane angular rotation [12, 13], we report in Sec. 
7 measurements with angular rotation of in-plane applied fields to test the theoretical 
predictions of the five models. We compare the experimental results with the predictions 
in Sec. 8, and we summarize our conclusions in Sec. 9. 

2. Generalized double-critical-state model (GDCSM) 

In an attempt to describe the fundamentals of the behaviour of a type-II superconductor 
subjected to time- varying fields and currents in situations where the current density J 
is not perpendicular to a theory was proposed in [14, 15, 16, 17, 18]. This approach 
has been referred to as the generalized double-critical-state model (GDCSM). 

Consider an infinite type-II superconducting slab centered on the yz plane. Suppose 
that in the absence of an applied current, an applied magnetic field produces a magnetic 
fiux density = Boao in the slab, where ao = ysinao + zcosaQ. Consider now what 
happens when the slab carries an average current density J. As a consequence of the 
self-field, vortices may be introduced, altering the magnetic fiux density inside the slab 
to B{x) — B{x)a^ where a{x) — y^ina{x) + zcosa(x). The current density J also 
depends upon x. 

Referring now to figure 1, let us resolve the current density and the electric field 
into components parallel and perpendicular io B{x)\ J = Jyy + J^z = J\\a -\- Jj_/3 
and E = Eyy + E^z = E\\a + £'±/3, where /3 = yco^a — z^ina = a x x. In high-/^ 
superconductors when is sufficiently above Sd, it is a good approximation to take 
B — iiqH inside the superconductor. In the steady state. Ampere's and Faraday's laws 
( J = V X H, V X = 0, with all quantities independent of time) then yield [18] 



J\\ — J cos (j) 
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Figure 1. Sketch of the unit vectors ^, i, a, and /3, the magnetic induction and 
the current density J and its components Jy and J± along a and /3. (/) is the angle 
between J and All vectors are parallel to the yz plane. 

^„^ + f^=0, (5) 

dec diXi 

The GDCSM was based on the fundamental ideas that a nonvanishing can 
occur only when the magnitude of is above Jc±5 the critical current density for flux 
depinning, and that a nonvanishing £"11 can occur only when the magnitude of J\\ is 
above Jc\\ , the critical current density for flux cutting. It was assumed in [18] that the 
superconducting material was isotropic and that the pinning centers were isotropically 
distributed, such that Jc± and Jc\\ could depend upon the temperature T but only upon 
the magnitude B of the magnetic induction B inside the sample. For thin samples, the 
current-induced self-field is much smaller than the applied magnetic induction Sq, so 
that Jc± and Jc\\ were taken as constants independent of x. The possible complications of 
surface pinning or surface barriers were ignored. One of the most important assumptions 
of the GDCSM was that fiux depinning and fiux cutting do not affect each other; i.e., it 
was implicitly assumed that Jc± was independent of J\\ and Jc\\ was independent of J^. 

Figure 2 provides a graphical display of the behaviour of the vortex array as a 
function of J±^ = J sine/) and Jy = J cose/) according to the GDCSM. In the middle 
of the open rectangle, the zone, for which < Jc± and |J||| < Jc\\^ neither fiux 
transport (depinning) nor fiux-line cutting occurs; i.e., both E± and £"11 are zero. Flux 
transport (without fiux-line cutting) occurs only in the zones labelled by T+ {Jj_ > Jc± 
and E± > 0) or T_ ( < —Jc± and E± < 0). Flux-line cutting (without fiux transport) 
occurs only in zones labelled by C+ ( Jy > Jc\\ and E\\ > 0) or C_ ( Jy < —Jc\\ and ^y < 0). 
Simultaneous fiux transport {\J±\ > Jc± and \Ej_\ > 0) and fiux-line cutting (|Jy| > Jc\\ 
and \E\\ \ > 0) occurs only in the zones labelled by C+T+, C_T+, C_T_, and C+T_. 




Figure 2. Behaviour of a vortex array of flux density B (which is parahel to the Jy 
axis) as a function of the magnitude J of the current density J and the angle ^ of J 
relative to B according to the generalized double-critical-state model [14, 15, 16, 17] 
(GDCSM), in which Jc± and Jc\\ are both assumed to be constants. 



According to the GDCSM, for given values of Jc± and Jc|h the critical current 
density Jc at the first onset of a nonvanishing electric field depends upon the angle 0. 
For = 7r/2 (J -L B)^ the critical current density for the onset of an electric field is 
the usual depinning critical current density, such that Jc = Jc±- On the other hand, 
for = (J II S), the critical current density for the onset of an electric field is the 
current density at the threshold for flux-line cutting, and Jc = Jc\\- For other values of 
0, we have 

Jc = T-^^ I tan 01 > tan (7) 
I sin0| 

= 1 — ^TT, I tan 01 < tan (/)c, (8) 
I cos 01 

where tan 0c = Jc±/Jc\\ • As is evident from figure 2 and (7) and (8), this model predicted 
that Jc as a function of has cusplike maxima at = ±0c and = ±(7r — 0c), where 



Jcmax = Y ^ii + Jc\\' example, the dashed curve in figure 3 maps the magnitude of 
J as a function of as the tip of the vector J traces the boundary of the zone shown in 
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Figure 3. Plot of Jc (in units of Jci_) vs (j) for two theoretical models accounting for 
flux-line cutting when Jc\\lJc^ = 2.5, for which (/)c = tan-i(0.4) = 0.381 rad = 21.8°. 
The dashed curve shows Jc/Jc± according to the generalized double critical-state model 
(GDCSM) [18], (7) and (8), and the solid curve shows Jc/Jc± according to the elliptic 
critical-state model [22, 23, 24], (11) and (12). 



figure 2. The rectangular shape of the zone in figure 2 arises from the assumptions of 
the GDCSM [14, 15, 16, 17, 18] that (a) the threshold Jc± for depinning is independent 
of the current density J\\ parallel to the vortices and (b) the threshold Jc\\ for fiux-line 
cutting is independent of the current density perpendicular to the vortices. 

Another prediction of the GDCSM [18] was that, in accord with figure 2, 
measurements of the angle of the electric field E just above the critical current as a 
function of the angle between J and Bq would show sharp changes at |0| = (/)c = 
tan~-^(Jc±/-^c||) and = n — (^c- If J is constrained by geometry to flow only in 
the z direction, then just above the critical current density, where J = JcZ and B = 
Bo(y sin ao + z cos ao)^ we have = —ao^ so that E = Eyy + E^z = E^^a + where 
a — ysinao+z cos ao and /3 = y cos ao — zsin ao. The GDCSM prediction was that when 
(j)c < ao < TV — (j)c 01 —7T + (j)c < ao < — the behaviour for J > Jc would be dominated 
by flux transport, such that E\\ — 0^ E 1. Bo^ E — E^jS — E_i{y cos ao — zsinao)^ and 

Ey/E, = -cotao. (9) 

On the other hand, when —(/)c < ao < (j)c^ the GDCSM prediction was that the 
behaviour for J > would be dominated by flux cutting, such that E±_ = 0, £^ || Sq, 
E — E\\a = E\\ [y sin ao + z cos ao)^ and 

Ey/Ez = tan^o- (10) 
The dashed curve in figure 4 shows the behaviour of Ey/Ez predicted by the GDCSM. 

3. Brandt and Mikitik model (BMM), extending the GDCSM model 

Brandt and Mikitik [19] proposed an extension (here called the BMM) of the GDCSM, 
in which the threshold for fiux-line cutting is reduced as a function of the magnitude of 
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Figure 4. Plots of the electric field ratio Ey/Ez vs ao for five theoretical models 
accounting for flux-line cutting. The predictions of the generalized double critical- 
state model and the Brandt-Mikitik model [GDCSM & BMM, dashed, (9) and (10)] 
and the original elliptic critical-state model [OECSM, dot-dashed, (20)] are all shown 
for tan(/)c = Jc^j Jc\\ = 0-4, for which 0c = tan~^ (Jc^/J^n) = 21.8°. Also shown is the 
behaviour of Ey/Ez predicted by the extended elliptic critical-state model [EECSM, 
solid, (24)] for an assumed value of r = 0.04, and that of the BLRM [dotted, (25)] for 
tan^ (/)e = (Je±/Jc||)' = 0.02. 



J^i and the threshold for depinning is reduced as a function of the magnitude of J\\. If 
these thresholds were bent enough, the shape of the surface of Jc(0) could be made to 
resemble an ellipse, as suggested by figure la of [19], thereby eliminating the cusplike 
behavior predicted by the GDCSM. In fact, if the BMM Jc(0) curve had the shape of 
an ellipse, it would be in exact agreement with the solid curve shown in figure 3. 

However, the BMM retained the assumption of the GDCSM that the electric field E 
is parallel to the J\\ axis on the flux-cutting-threshold portions of the Jc(0) curve closest 
to the J|| axis and that E is parallel to the axis on the flux-transport-threshold 
portions of the Jc(0) curve closest to the axis. Thus, in common with the GDCSM, 
the BMM predicts a 90 degree change in the direction of E when (f) passes through the 
point of intersection of the curves representing the thresholds for flux cutting and flux 
transport. If this intersection occurs when \ J^lJ\\\ — Jc^lJcW-) then the shape of the 
curve of Ey/Ez curve agrees exactly with that of the GDCSM, as shown in figure 4. 

4. Original elliptic critical-state model (OECSM) 

Experimental results of Fisher et al [20, 21] on the time evolution of the static magnetic 
moment of a superconducting plate subjected to an alternating magnetic field applied 
perpendicular to a dc magnetic field were found not to be well described theoretically 
by calculations based on the GDCSM. However, Romero-Salazar and Perez-Rodriguez 
[22, 23, 24] introduced an elliptic critical-state model, which they found yielded good 
theoretical agreement with the experiments of [20, 21]. 

Figure 5 provides a graphical display of the behaviour of the vortex array as a 
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function of = J sine/) and Jy = J cose/) according to the elliptic critical-state model 
[22, 23, 24]. We interpret this behaviour as follows: In the zone inside the ellipse 
described by 

sin^ (/) cos^ 1 



J,^{Bf J,||(5)2 Je(5,0)2 



or 



sin^ 6 cos^ ' 



= + (12) 

neither flux transport (depinning) nor flux cutting occurs {E^ — and E\\ — 0). Flux 
transport, for which the vortices are depinned {E^ ^ 0), occurs everywhere outside the 
ellipse (except when = 0) in zones with labels including the symbol {E± > 0) or 
T_ {E± < 0). Flux-line cutting, for which £"11 7^ 0, occurs everywhere outside the ellipse 
(except when J|| = 0) in zones with labels including the symbol C+ {E\\ > 0) or C_ 
{E\\ < 0). In this model, both flux transport and flux-line cutting occur simultaneously 
nearly everywhere outside the ellipse of Jc(0) vs (/). The solid curve in flgure 3 gives an 
example of the behaviour of Jc vs (j) using the elliptic critical-state model. 

We note here that the behaviour at the critical current according to the elliptic 
critical-state model (see figure 5) requires that fiux depinning and fiux cutting be 
interdependent. To formulate this, we introduce the following new notation for the 
critical current densities at the thresholds for flux depinning Jcd{B^ (j)) and flux cutting 
Jcc{B^(/))^ where (/) is the angle between J and as in flgure 1. Our interpretation 
of the underlying physics of the elliptic critical-state model [22, 23, 24] is that (a) the 
threshold Jed for depinning is a monotonically decreasing function of the magnitude of 
J||, decreasing from Jc± when Jy = to zero when |J||| = Jc\\^ and (b) the threshold 
Jcc for flux-line cutting is a monotonically decreasing function of the magnitude of J^, 
decreasing from Jc\\ when = to zero when | = Jc±- In summary, at the critical 
current density, J = Jc{B^ 0), and J lies on the ellipse of (12) as shown in flgure 5, such 
that |J^| = Jcd{B,T) = J,(S,0)|sin(/)|, |J||| = Jcc{B,T) = J,(S, 0)1 cos(/)|, and 

V sm + tan^ 0c cos^ 
where tan 0c = Jc±{B) / Jc\\{B) . Note that within this model Jc±{B) = Jc{B,7r/2) and 

Jc||(5) = Jc(5,0). 

Further fundamental assumptions of the original elliptic critical-state model 
[22, 23, 24] are that the electric fleld obeys 

E(J) = 0, J<Jc{B,4>), (14) 

= p[J-JciB,4>)], J>UB,4>), (15) 

and that the components of E perpendicular and parallel to B obey [22, 23, 24] 

E±/E = J^/Jex, (16) 
E^l/E = J||/Jc||. (17) 




Figure 5. Behaviour of a vortex array of flux density S as a function of the magnitude 
J and angle (j) of the current density J relative to the direction of B when flux-line 
cutting and depinning interact and the critical current density is given by the elliptic 
critical-state model of [22, 23, 24]. Here J is shown at the critical current density, 
where J = Jc^ Jcc = Jc coscj) is the threshold for flux cutting, and Jed = Jc sin(/) is the 
threshold for flux transport (depinning). 



If the current is constrained by geometry to flow only in the z direction, then at the 
critical current density, when J = zJc and B = Bo{y sin ao+z cos q^o), we have = — o^o, 
such that J|| = JcCosao and J± = — JcSinao (see figure 1) and 

E^/E= -(J,/J,^)sinao, (18) 

Eii/E = (J,/J,||)cosao. (19) 
Since E = Eyy + E^z = E\\a + Ej_/3 and E\\/Ej_ = —{Jc±/ Jc\\)/ tano^o, we have 

Ey _ {E\\/E^)t^nao + 1 _ (tan^^ - l)tanao 

Ez (£'||/£'^) — tano^o tan (/)c + tan^ o^o 
where tan^c = Jc^/Jc\\' Experiments generally yield tan^c < 1. According to 
this form of the elliptic critical-state model, as a function of tano^o, Ey/E^ has 
extrema (a maximum and a minimum) when tano^o =FA/taii0^, where Ey/E^ 
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±|1 — tan (j)c\/2^/i^^n^c' The dot-dashed curve in figure 4 shows the behaviour of Ey/E^ 
predicted by the original eUiptic critical-state model [22, 23, 24], (20). 

5. Extended elliptic critical-state model (EECSM) 

However, the original elliptic critical-state model [22, 23, 24] suffers from an important 
deficiency. By summing the squares of (16) and (17), we obtain 

rj rj Jci_ Jc\\ 

The right-hand side can be equal to one only when J — Jc^ SiS can be seen by making 
use of (18), (19), and (11). However, when J — Jc^ we must have E — 0. Therefore, 
additional equations are required in order to apply the elliptic critical-state model to 
finite values of the electric field when J > Jc^ where the dynamical processes of fiux 
transport and fiux cutting become important. 

We therefore propose an extension of the elliptic critical-state model to account for 
dissipative processes when J > Jc hy writing 

Ei_ = Jx, (22) 

^11 -PcJ\\. (23) 

where p/ is a nonlinear function of J, S, and T with a positive value describing 
fiux fiow and pc is another nonlinear function of J, S, q^q? and T with a positive value 
describing fiux cutting. In the absence of fiux creep, is zero when | is less than 
the threshold Jed for depinning, and Pc is zero when |J||| is less than the threshold Jcc 
for fiux cutting. To clarify what we mean by a threshold for fiux-line cutting, consider 
starting with the sample carrying a current for which Jy = and < Jj_ < Jc_l, such 
that E — 0. If J|| is increased above the threshold value Jcc, fiux-line cutting begins 
{E\\ > 0), but the elliptic critical-state model (and the experiments to be presented 
later) show that this is accompanied by fiux transport {E_\_ > 0). In other words, when 
J > Jc{(l))^ fiux cutting and fiux transport occur simultaneously; the thresholds for fiux 
cutting and fiux transport are intimately linked. 

If the current is constrained by geometry to fiow only in the z direction, then 
at the critical current density, when J = zJ and B = Bo^ys'mao + zcosao), we 
have (j) = — Q^o, such that J\\ = Jcoso^o and J± = — Jsino^o (see figure 1). Since 
E = Eyy + EzZ = E\\a + Ej_/3^ where a — y8'mao + z cos and /3 = y cos ao — z sin ao^ 
and E\\/E±_ = —{Pc/Pf)/ tano^o^ we obtain, making use of (22) and (23), 

Ey_ _ {E\\/E^)imiaQ + 1 _ (r - 1) tanap 
Ez {E\\/E^) — tan r + tan^ 

where now r = pd Pf^ which must be a positive number. 

Note that (24) has the same mathematical form as (20), except that tan 0c is 
replaced by r. According to this extension of the theory, EyjE^ does not depend 
at all upon tan 0c — Jc^/Jc\\ but instead depends upon r, the ratio of the effective 
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flux-cutting resistivity to the effective flux-flow resistivity. If r is independent of 
tano^o, Ey/Ez has extrema (a maximum and a minimum) when tano^o — TV^? where 
Ey/Ez = ±|1 — '^l/Sy^- The parameter r thus plays a crucial role in describing 
the dynamical properties of the mixed state above the critical current density Jc- 
Unfortunately, experiments to date typically do not report values of pc and p/, from 
which r = Pc/pf could be independently measured. 

The solid curve in flgure 4 shows the behaviour of Ey/E^ predicted by this extended 
elliptic critical-state model (24) for an assumed value of r = 0.04. 

6. Badia-Majos, Lopez, and Ruiz model (BLRM) 

In [25], Badia-Majos et a/ argued that for J just above Jc{(p)^ the direction of the electric 
field should be perpendicular to the Jc(0) curve. Under the same conditions leading to 
Eqs. (20) and (24), their model (here called the BLRM) leads to the expression 

Ey 1 dJcjap) (tan^ (/)c l)tanaQ 

Ez Jc{<^o) dao tan^ (/)c + tan^ 

where the latter equation applies when Jc(0) is given by the ellipse of (12) (see figure 
5) and tan^c = Jc±/Jc\\- The dotted curve in figure 4 shows the behaviour of Ey/E^ 
predicted by the BLRM for an assumed value of tan^ = {Jc±/Jc\\Y — 0.02. 

7. Experiment 

We next report experiments to test the fundamental assumptions of the above-described 
theories. While the predictions have been made for an infinite superconducting slab, 
we have carried out experiments on a current-carrying superconducting film of finite 
dimensions subjected to an in-plane magnetic field. We have assumed that the effects of 
finite dimensions and self-fields are small enough that the theoretically predicted effects 
dominate and the experiments can determine which of the proposed theories is correct. 

A sketch of the measured device, a YBCO film of thickness 500 nm, is shown in 
figure 6. The track is of width Wy = 200 /xm, and two sets of voltage taps are positioned 
along its length, = 3000 pm apart from each other. The voltage/current taps are 
labelled A, B, C, D, E, F, G, H, K, L, A', B'. A two- axis goniometer [26] was employed 
to measure the longitudinal and transverse voltage simultaneously by a four-terminal 
technique, using the two channels of the employed voltmeter, while a magnetic field was 
swept in-plane. For all measurements the current / was injected through the taps A 
and B and taken out through A' and B', such that the current direction was (A, B) 
(A', B'), i.e., in the positive z direction. The current was applied in short pulses 
to minimize sample heating. Longitudinal and transverse voltages were measured using 
combinations of the taps C, . . . , L. 

The applied field was oriented to be in the plane of the film, such that angle 9 
between the applied field and the film normal ft was as close as possible to 90°, as 
described below. The current density J was constrained to flow in the z direction. 
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Figure 6. The pattern used. The track is shown in black and the volt age /current 
taps in gray with their labels A, B, C, D, E, F, G, H, K, L, A', and B'. Voltage 
measurements were made over the length Wz = 3000 /im and width Wy = 200 jam. 

while the angle a = ao of the applied magnetic induction Bq was oriented as shown in 
figure 1, such that the angle of J relative to Bq was = — o^o- 

In order to align the sample such that Bo\\ab^ measurements of Jc{0) were performed 
at = ±90°, corresponding to maximum Lorentz force {Bq _L J). The angle 9 = 90° 
(where Bo\\ab) was assumed to be reached where Jc{0) reached a maximum, and the 
following in-plane scans were performed with this value of 9. 

For the reported measurements, Jc was determined from Vz{I) curves only, where 
Vz is the voltage measured along the length of the track (e.g., by using taps C and G). 
For each value of ao, Ic was determined as the current where reaches the selected 
voltage criterion Vc. The transverse voltage Vy{I) was measured simultaneously (e.g., 
using taps G and H). The voltages were then converted to electric fields [E^ = Vc/wz 
and Ey = Vy{Ic)/wy] in order to obtain Ey/Ez as a function of aQ. 

Figure 7 presents plots of Ey/Ez vs ao obtained at T = 85 K and /hqH = 1 T 
for voltage criteria K = 15 /xV (square symbols), 30 /xV (dots), 60 /iV (triangles), and 
150 /xV (inverted triangles). Note that the scatter in the values of Ey/Ez becomes 
smaller with increasing values of V^. In the raw data, Ey/Ez reached zero at a value 
of Q^o = +3°. This could be understood by the fact that the alignment of the in-plane 
angle was performed by visual inspection prior to sample mounting. Consequently, we 
shifted the curves by —3° in q^q, so that Ey/Ez = at o^o = 0° in all the plots shown 
here and later. The solid curve in figure 7 shows the function Ey/Ez = — coto^o (9), 




Figure 7. The measured ratio Ey/Ez vs ao at T = 85 K and ji^H = 1 T. The 
sohd curve is Ey/E^ = — cotofo (9), expected due to flux transport in the absence of 
flux-hne cutting. 




Figure 8. The measured ratio Ey/E^ vs at T = 77.35 K and iiqH = 8 T. The 
sohd curve is Ey/E^ = — coto^o (9), expected due to flux transport in the absence of 
flux-line cutting. 



the result that would be expected for no flux-line cutting; i.e., when £"11 = 0, _L Sq, 
and E = Eyy + E^z = E±/3 = E±{y cos ao — zsino^o). (Note from figure 1 that in this 
experiment (/) = — o^o, so that when ao > 0, we have E^ > 0^ Ey < 0, and Ej_ < 0.) 

Similar measurements were also performed at T = 77.35 K and /hqH = 8 T for 
Vc = 15 /xV (square symbols) and 60 /iV (triangles), and the results are shown in 
figure 8. Note that, as in figure 7, the scatter in the values of Ey/E^ becomes smaller 
with increasing values of the voltage criterion Vc. 

The data presented in Figs. 7 and 8 show results measured at constant Ez^ but 
we also did an analysis making use of the saved V{I) data of both voltmeter channels 
to construct corresponding plots of Ey/Ez vs at constant Etot = ^^i^^ 
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Figure 9. Critical-current density Jc (K = 60 /iV, Ez = 20 mV/m) vs the angle ao 
of the in-plane applied field at 85 K and 1 T (inverted triangles) and 77.35 K and 8 
T (open circles). The corresponding solid curves are fits using the elliptic model (12) 
with Jell = 2.65 X 10^ A/m^ and Jc± = 0.895 x 10^ A/m^ at 85 K and J^y = 7.86 x 10^ 
A/m^ and Jc± = 2.06 x 10^ A/m^ at 77.35 K. 

there were some differences from the results for constant for the smallest values of 
these were not significant for the larger values of K- 

8. Comparison with theoretical predictions 

Shown in figure 9 are the critical current densities Jc vs ao measured at T = 85 K and 
IHqH = 1 T (inverted triangles) and at T = 77.35 K and /hqH = 8 T (open circles), 
both with the voltage criterion Vc = 60 /xV, which corresponds to = 20 mV/m. The 
corresponding solid curves are fits to the data using the elliptic critical-state model: at 
T = 85 K and fioH = 1 T, J^y = 2.65 x 10^ A/m^ and Jc± = 0.895 x 10^ A/m^ (the 
average of values at ±90°), such that Jc±/Jc\\ = 0.338 ± 0.002, and at T = 77.35 K and 
IUqH = 8 T, Jell = 7.86 X 10^ A/m^ and Jc± = 2.06 x 10^ A/m^ (the average of values at 
±90°), such that Jc±/Jc\\ = 0.262 ± 0.004. (The indicated errors are chosen such that 
both values lie within the range.) The data are rather well fitted by the elliptic critical- 
state model as in [27]. In agreement with previous experiments [28, 29, 30, 31] the 
data clearly do not show the cusp-like maxima predicted by the GDCSM (see figure 3). 
Thus the Jc vs ao measurements support either the original elliptical critical-state model 
(OECSM) or the extended elliptical critical-state model (EECSM), since both assume 
the same dependence of Jc{(fy) but do not support the generalized double critical-state 
model (GDCSM). 

Figure 10 shows comparisons of the measurements of Ey/E^ vs at 85 K and 1 
T (using the voltage criterion = 60 /xV) with the predictions of the original elliptic 
critical-state model [dot-dashed curve, (20)] and the Badia-Majos-Lopez-Ruiz model 
[dotted curve, (25)] using the measured tan^c = Jc^j Jc\\ = 0.338 , the extended elliptic 
critical-state model using r = pd Pf = 0.090 assumed to be independent of [solid 
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Figure 10. Comparison of the measured Ey/Ez (open symbols) vs ao at 85 K and 1 T 
{Vc = 60 //V) with theoretical curves for the OECSM [(20), dot-dashed] and the BLRM 
[dotted, (25)] with tan 0c = Jc±/Jc\\ = 0.338, the EECSM with r = pc/pf = 0.090 
[(24), solid], and the GDCSM and BMM with tan 0c = Jc±/Jc\\ = 0.338 [(9) and (10), 
dashed] . 

curve, (24)], and the generalized double critical-state model and the Brandt-Mikitik 
model using the measured tan^c = Jc±/Jc\\ = 0.338 [dashed curve, (9) and (10)]. The 
theoretical curves for the extended elliptic critical-state model and the Badia-Majos- 
Lopez-Ruiz model can account for the large values for the magnitude of the measured 
Ey/Ez, while the original elliptic critical-state model cannot. The data also clearly do 
not show the sawtoothlike features predicted by the generalized double critical-state 
model and the Brandt-Mikitik model. To summarize, these measurements of Ey/Ez 
vs Q^o favor the extended elliptical critical-state model and the Badia-Majos-Lopez-Ruiz 
model but do not support the original elliptic critical-state model, the generalized double 
critical-state model, or the Brandt-Mikitik model. 

Figure 11 shows comparisons of the measurements of Ey/Ez vs ao at T = 77.35 K 
and 8 T (using the voltage criterion Vc = 60 /xV) with the predictions of the original 
elliptic critical-state model [dot-dashed curve, (20)] and the Badia-Majos-Lopez-Ruiz 
model [dotted curve, (25)] using the measured tan^c = Jc^/ Jc\\ — 0.262, the extended 
elliptic critical-state model using r = Pc/Pf = 0.026 assumed to be independent of 
Q^o [solid curve, (24)], and the generalized double critical-state model and the Brandt- 
Mikitik model using the measured tan^c = Jc^/Jc\\ = 0.262 [dashed curve, (9) and 
(10)]. The theoretical curve for the extended elliptic critical-state model can partially 
account for the large values for the magnitude of the measured Ey/Ez, but the original 
elliptic critical-state model and the Badia-Majos-Lopez-Ruiz model cannot. The large 
magnitudes of the experimental values of Ey/Ez at 77.35 K and 8 T cannot be 
understood with any of the proposed models without modification. 

Although we cannot at present provide a justification for the following assumption. 
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Figure 11. Comparison of the measured Ey/Ez (open symbols) vs at 77.35 K and 
8 T (Vc = 60 iiY) with theoretical curves for the OECSM [(20), dot-dashed] and the 
BLRM [dotted, (25)] with tanc/)^ = Jc^/Jc\\ = 0.262, the EECSM with r = pc/pf = 
0.026 [(24), solid], and the GDCSM and BMM with tan 0c = Jc±/Jc\\ = 0.262 [(9) and 
(10), dashed]. 
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Figure 12. The measured Ey/E^ (open symbols) vs ao at 77.35 K and 8 T 
(Yc = 60 /iV, Ez = 0.02 V/m) and the assumed fit function [(26), solid curve] with 
the fit parameters r = Pc/pf = 0.080 and c = 1.936) 

if we assume the following form for theoretical values of Ey/Ez, 
E, ^ ^(r-l)tanao 
Ez r + tan^Q^o 

we obtain a considerably better fit to the data at 77.35 K and 8 T, as shown in figure 
12, where the fit parameters were found to be r = 0.080 and c = 1.936. 
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9. Conclusions 

The purpose of this paper has been to compare five theoretical models describing the 
critical state of type-II superconductors when the current density J is at an arbitrary 
angle with respect to the local magnetic flux density B and then to use experimental 
measurements of the angular dependence of the critical current density Jc and the 
electric fleld E (for J just above Jc) as a test of these theories. 

For simplicity, the geometry of an infinite, macroscopically isotropic type-II 
superconducting slab of finite thickness centered on the yz plane was used for calculating 
the predictions of the five theories. However, the experiments were performed using a 
thin YBCO film subjected to in-plane currents and magnetic fields, but we assumed that 
the effects of finite dimensions and self-fields were small enough that the experiments 
would be able to make a clear distinction among the five theories. 

Measurements of the angular dependence of the critical-current density Jc showed 
that the cusp-like behaviour predicted by the generalized double critical-state model 
(GDCSM) was not verified. Instead, the angular dependence of Jc showed a behaviour 
very similar to that assumed by both the original elliptical critical-state model (OECSM) 
and the extended elliptical critical-state model (EECSM). 

Measurements of the angular dependence of the ratio of the transverse to the 
longitudinal components of the electric field Ey/ at 85 K and 1 T for J just above 
Jc also showed that the sawtoothlike behaviour predicted by the generalized double 
critical-state model (GDCSM) and the Brandt-Mikitik model (EMM) was not verified. 
While all three elliptical critical-state models predict a smooth angular dependence 
of Ey/Ez^ as found experimentally, the magnitude of Ey/E^ predicted by the original 
elliptical critical-state model (OECSM) was too small to account for the experimental 
values. On the other hand, both the extended elliptical critical-state model (EECSM) 
with r = Pc/pf = 0.090 and the Badia-Majos-Lopez-Ruiz model (BLRM) were found to 
be in reasonable agreement with the data at 85 K and 1 T. However, the EECSM has 
the advantage that it is applicable for currents J well above Jc{(fy'> while the BLRM is 
currently limited to J just above Jc{(fy' 

None of the five models was found to provide a good description of the data for 
Ey/Ez at 77.35 K and 8 T. A small value of r = pc/ Pf = 0.026 within the extended 
elliptical critical-state model (EECSM) can account roughly for the large magnitudes 
of Ey/Ez^ but the resulting theoretical curve does not provide a good fit to the shape 
of the data. By adding an additional fit parameter to the theoretical expression for 
Ey/Ez within the EECSM, better agreement with the magnitude and shape of the 
data for Ey/Ez at 77.35 K and 8 T can be obtained, but at present we have no 
theoretical justification for this assumption. Further theoretical and experimental work 
to investigate the influence of flnite dimensions and self-flelds will be needed to shed 
light on this problem. Our theory assumes homogeneity of the electric fleld but this 
needs to be conflrmed, since earlier experiments measuring the local electric fleld in 
the longitudinal geometry found strong inhomogeneities in E along the sample length 
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[33, 34, 35, 36, 37]. 

The experimental results in this paper call attention to two important features 
that to date have not been generally recognized. First, when J is neither parallel 
nor perpendicular to the local magnetic flux density both flux cutting and flux 
transport occur simultaneously when J exceeds the critical current density This 
indicates the intimate relationship between flux cutting and depinning. Second, the 
dynamical properties of the superconductor when J exceeds Jc depend in detail upon 
two effective resistivities, pc and p/, and their ratio r = Pc/ Pf- The consequence of this 
is that to provide a full description of the properties of type-II superconductors in the 
resistive state, both p^ and pf should be measured. It can be shown that in cylindrical 
geometry, values of r < 1 are responsible for the paramagnetic moment generated along 
the length of a current-carrying type-II superconductor above its critical current when 
it is subjected to a longitudinal applied magnetic held [8]. 

It seems likely that the angular dependence of the critical current density in type-II 
superconductors will depend upon the details of the underlying pinning structure, and it 
is possible that in many cases Jc{(p)^ which can be obtained experimentally, will not be 
given precisely by the elliptic form given in (11) and (12). For example, Jc(0) niay even 
take the form of a superellipse, as suggested in [25] and [32]. In other words, although 
one should not take the elliptic model too seriously, it seems likely that Jc(0) generally 
has roughly the shape shown in flgure 5, and that it is a general principle that both flux 
cutting and flux transport occur simultaneously when J > Jc{(t))' 
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